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' We study radial solutions in a ball of M™ of a semilinear, parabolic-elliptic 

Patlak-Keller-Segel system with a nonlinear sensitivity involving a critical power. 
For N = 2, the latter reduces to the classical "linear" model, well-known for its 
critical mass Stt. We show that a critical mass phenomenon also occurs for 

■ Af > 3, but with a strongly different qualitative behaviour. More precisely, if 

■ the total mass of cells is smaller or equal to the critical mass M, then the cell 
' density converges to a regular steady state with support strictly inside the ball 

as time goes to infinity. In the case of the critical mass, this result is nontrivial 
since there exists a continuum of stationary solutions and is moreover in sharp 

■ contrast with the case JV = 2 where infinite time blow-up occurs. If the total 

mass of cells is larger than M, then all solutions blow up in finite time. This 
actually follows from the existence (unlike for A'^ = 2) of a family of self-similar, 
blowing up solutions with support strictly inside the ball. 
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1 Introduction 

1.1 Origin of the problem 

Chemotaxis is the biological phenomenon whereby some cells or bacteria direct 
their movement according to some chemical present in their environment which 
can be attractive or repulsive. We shall focus on the case where the chemical is 
attractive (tlien called chemoattractant) and self-emitted by cells. For instance, 
in case of starvation, amoebas Dyctyostelium discoideum emit cyclic adenosine 
monophosphate (cAMP) which attract themselves. Chemotaxis is thus a strong 
mean of communication for cells and leads to collective motion. 
For more details on the social life of amoebas Dyctyostelium discoideum, see the 
article [17] of M.A. Herrero and L. Sastre. 



1.1.1 Mathematical formulation 

Assuming that cells and chemoattractant are diffusing and that cells are sensi- 
tive to the chemical's concentration gradient (a fact experimentally observed), 
Patlak in 1953 (cf. [35]) and Keller and Segel in 1970 (cf. 24|) have pro- 
posed the following mathematical model, a parabolic-parabolic system known 
as Patlak-Keller-Segel system : 

Pt = D,Ap^V[xVc] (1) 
Ct^D2Ac + fip (2) 

where p is the cell density, c the chemoattractant concentration, Di and D2 are 
diffusion coefficients, x = x(P: c) is the sensitivity of cells to the chemoattractant 
and n the creation rate of chemical by cells. 

Cells and chemoattractant are assumed to lie in a bounded domain il of 
with N >2 (M'^ or M.^ physically speaking) so we have to specify the boundary 
conditions. 

For the cell density p, it is natural to impose a no flux boundary condition 

D^^-X{p,c)^=0 ondn, (3) 
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where v denotes the outward unit normal vector to the boundary d^. 

For the chemoattractant concentration c, a Dirichlet or Neumann boundary 

condition is assumed : 

3c 

c = on af7 or = on d^. (4) 

ov 

In the case of "fat" cells, the chemoattractant diffuses much faster than cells. 
Hence, two timescales appear in the system and to the limit, we can assume that 
the chemical concentration c reaches instantaneously its stationary state. After 
renormalization, these considerations lead to the simplified parabolic-elliptic 
Patlak-Keller-Segel system : 

Pt = Ap-V[xVc] (5) 
-Ac = p (6) 

with the same boundary conditions as above, which then become : 

|^-X(p,c)|^=0 onaf7. (7) 
3c 

c = on 9f7 or — = on d^. (8) 

ov 

We would like to add that this hypothesis of " fat" cells seems all the more ap- 
propriate than some smaller cells like Escherichia coli have a 'run and tumble' 
motion for which diffusion does not seem to be the most suitable modeling. On 
that subject, see the article of B. Perthame [32] for a kinetic approach which 
takes into account these characteristics and allows to recover the Patlak-Keller- 
Segel model in a diffusion limit. 

For a review on mathematics of chemotaxis, see the chapter written by M.A. 
Herrero in [TB]. For a review on the Patlak-Keller-Segel model, see both articles 
of D. Horstmann [TH] and pU] . 

In [21], D. Horstmann and M. Winkler have studied the case where the sensi- 
tivity X depends only on p and shown that : 

_2 

N 

and is even uniformly bounded in time 



• If x(p) < C'p' for p > 1 and q < then the cell density p exists globally 



• If x(p) > Cp^ for p> 1 and q > jj, then p can blow up. 
See also [H US] [Ml 1201 13S] for related results. 

Thus, the power q = jj oi the nonlinearity x{p) is critical for that system. This 
is why we are going to focus on the following problem, noted [PKSq), with a 
special interest to the case q — : 

Pt = Ap-Vb^Vc] (9) 
-Ac^p (10) 

where the boundary conditions become 

1^ - 1^ - on dn. (11) 

ov ov 

c = ondn or = on d^. (12) 

ov 
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Remark 1.1 System (PKSq) can be seen as the macroscopic description of 
a collection of n particles following a generalized stochastic Langevin equa- 
tion. Making a mean field approximation, it is actually obtained as a nonlinear 
Fokker-Planck equation in a proper thermodynamic limit n — >■ cx). For more 
details, see the article of P.H. Chavanis fl(A section 3.4-]. 



1.1.2 Radial setting 

In this paper, we restrict our study to tlie case of radially symmetric solutions 
of {PKSq) where VI is the open unit ball B C centered at the origin. Note 
that by using the scaling of the system and its invariance by translation, we can 
of course cover the case of any open ball of R^. 

We would like to point out that for = 2, the critical exponent is q = 1, so the 
latter system reduces to the most studied Keller-Segel parabolic-elliptic type 
system : 

Pt^^p-W[pWc] (13) 
-Ac = p (14) 

It is a well-known fact that this system exhibits a critical mass phenomenon. 
More precisely, denoting m the total mass of the cells, it has been shown for 
radially symmetric solutions that : 

• If m < Stt, then p{t) is global and converges to a steady state as t goes to 
infinity. 

(sec 3 by P. Biler, G. Karch, P. Laurengot and T. Nadzieja). 

• If m = Stt, then p{t) blows up in infinite time to a Dirac mass centered at 
the origin. 

(see again [3] and [53] by N.I. Kavallaris and P. Souplet for refined asymp- 
totics). 

• If m > Stt, then p{t) blows up in finite time to a Dirac mass, 
(see [IH] by M.A. Herrero and J.L. Velazquez). 

Moreover, this system exhibits a similar phenomenon in the case of the whole 
plane R^. See the work of [H [SJ [71 [T3]. In the nonradial case in a bounded 
domain, results are slightly different (see the book 37 of T. Suzuki). The be- 
haviour of the parabolic-parabolic system in seems more intricate. See [Tllllj. 



From now on, we consider the case iV > 3. 



Adapting the procedure described in the article jSj of P. Biler, D. Hilhorst and 
T. Nadzieja (or also in [21 123]): we can reduce the system {PKS)q to a single 
one-dimensional equation. 

Indeed, denoting Q{t,r) the total mass of the cells in B{0,r) at time t for 
< r < 1, we formally get : 
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Qt Qrr 



N- 1 



Qr 



Qr 



CTN r 



N-l 



(15) 



where an denotes the "area" of the unit sphere in M^. 
Then, setting P{t,x) = Q{t,xT^), we obtain : 



Pt = N^X^-- Pr 



' N ■ 



(16) 
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Finally, setting u{t,x) — — P{-^,x) where V^r = ^ is the "volume" of B, 

Ni Vn 

we get : 

Ut^x'^~^U^a:+u{u^y. (17) 

Moreover, by the no flux boundary condition on the cell density, it is formally 
clear that the total cell mass m is constant in time. Hence, setting 



N^Vn 

we also have the boundary conditions that for all ^ > 0, 

u{t,0) = 

u{t, 1) = m. 
A simple calculation also shows that 

p{t,x) u^{NH,x). 

Hence, p is simply proportional to u^, up to a time rescaling. It means that 
the derivative of u is the quantity with physical meaning and should then be 
nonnegative. 

Eventually, we shall focus on the following problem, noted {PDE„i), which we 
consider as our model for chemotaxis : 

ut = x'^^ ^ Uxx + u {u^y t>0 0<a;<l (18) 

u{t,0) = t>0 (19) 

u{t,l)=m, t>0 (20) 

Uxit,x)>0 t>0 0<x<l. (21) 

Equation ([18)) presents two difficulties since the diffusion is degenerate at a; = 
and the nonlinearity is not Lipschitz continuous. We shall assume that the 
initial data uq belongs to the class 

Y,n — {u e C([0; 1]), u nondecreasing , u{0) exists, u{0) = 0, u{l) = m}. 

For such Uq, we have established in [29j the existence and uniqueness of a max- 
imal classical solution u such that u{t) S Ym for all t G [0,Tmaxiuo)), where 
Tmaxiuo) is the maximal existence time. See Subsection 12.11 below for precise 
definitions and more details. 

1.2 Main results 

We now focus on the case of the critical exponent q — jj with > 3. 
The set of stationary solutions can be precisely described. 

We shall prove that the stationary solutions of [PDEm] are the restrictions to 
[0, 1] of a family of functions {Ua)a>vi with the following properties : 

• C/i(0) = 0, Ui is nondecreasing and reaches its maximum M aX x = A 
from which Ui is flat. 

• All {Ua)a>Q are obtained by dilation of Ui, i.e. Ua{x) — Ui{ax) for all 
x>0. 
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Using this, we can then prove : 

Theorem 1.1 Let m > 0. Considering problem {PDE„i) with q = jj : 

i) If < m < M, then there exists a unique stationary solution. 

a) If m = M, there exists an continuum of steady states : (C^a|[o,i])(j>^- 
Note that the corresponding cell densities have their support strictly inside 
B. 

Hi) If m> M, there is no stationary solution. 
The previous theorem leads us to set the following definition. 
Definition 1.1 

We call M the critical mass of problem (PDEm) for q= jf and 

M = N^Vn X M 
the critical mass of system {PKSq). 

We would like to stress that the system (PKSq) exhibits two levels of criticality. 
We have already seen the first level which consist in choosing the right exponent 
q= in order to balance the diffusion and aggregation forces in the system. 
Another observation is that the norm of the cell density p is invariant by 
the scaling of the system exactly ioi q = j^. This fact opens the door to the 
possibility of a critical mass phenomenon. The following two theorems state 
that it indeed occurs. 

Theorem 1.2 Let N > 3 and q = . 

If m < M and uq G Y^, then 

u{t) Ua in C\[0,1]) 

t— >+CX) 

for some a > 0. 

Theorem 1.3 Let N > 3 and q = jf . 
Ifm> M, then for all uq S Y^, 

Tmax{uo) < 00. 

Moreover, denoting J\f[u] = sup we have : 

xeio,i] 

lim Af[u{t)] = 00. 

In addition, for slightly supercritical mass, we can show the existence of blowing- 
up self-similar solutions. 

Theorem 1.4 There exists M+ > M such that for all m € (M, M+] there 
exists a family of blowing-up self-similar solutions of problem (PDEm)- 
Moreover, the corresponding cell densities have their support strictly inside B. 
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1.3 Comments and related results 



1.3.1 Description of the ideas of the proofs 

The global existence part of Theorem ll.2l for subcritical or critical mass is based 
on comparison with suitable supersolutions, combined with some continuation 
results obtained in [53]. Our convergence statements heavily rely on Lyapunov 
functional type arguments. 

More precisely, we show that the evolution problem {PDEm) induces a gradient 
type dynamical system on = F„ n C^([0, 1]), with global relatively compact 
trajectories. Moreover, we exhibit a strict Lyapunov functional : 

J-\u\~ I — ax. 

Indeed, formally, it is easy to check that 



—T[u{t)\ = / Mf^-j ^dx = Ut- - / 

it Jo l-q x^ 1 L 1 - 'J'Jo -/o 

inditions and to (ITS)) , we t 

J-[u(<)] = - / .^""'f dx. 
^ ^ " Jo uix^-i 



Ut 



dx 1 — q 



dx. 



Thanks to the boundary conditions and to (ITS)) , we then have 

d_ 
dt" 

However, this computation is not rigorously valid, since can vanish on a 
whole interval for instance. Nevertheless, we can overcome this difficulty and 
prove that is indeed a strict Lyapunov functional by expressing it as the limit 
as e goes to zero of a family of strict Lyapunov functionals J^e for suitable ap- 
proximate problems (cf. problem {PDE!^^^) introduced in Subsection l2.2.2p . We 
note that the proof of the compactness of the trajectories relies on a different 
transformation, leading to another auxiliary problem (tPDEm) (cf. Subsection 
12.2.11 below). In the subcritical case, since there is a single steady state, this 
immediately implies the convergence of the trajectory. But in the critical case, 
the situation is more delicate, since there exists a continuum of steady states 
and the solution could oscillate without stabilizing. However, thanks to a good 
relation between order and topology of the set of stationary solutions, we can 
prove stabilization by arguments in the spirit of (though simpler than) those in 
the articles ED [35] of H. Matano and T.I. Zelenyak. 



As for our blow-up results (Theorems 11.31 and II. 4p . their proofs relies on the 
construction of a subsolution which becomes a self-similar solution after some 
time. The latter's profile is solution of an appropriate auxiliary ordinary dif- 
ferential equation which is a perturbation of the stationary solution's equation. 
The construction, as well as the study of the steady states (cf. Theorem II. IL 
requires some rather delicate ODE arguments. 



1.3.2 Open problems 

i) A natural and very interesting question would be to determine the basin 
of attraction of a given steady state Ua with a> A. 

ii) For the self-similar solutions in Theorem ll.4l it is easy to see that the blow- 
up rate of the central density of cells (proportional to Ux{t, 0)) behaves like 
(Tmax — ■ It would be interesting to know if all solutions of problem 
{PDEm) blow up at the self-similar rate or if there also exists blow-up of 
type //. 
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1.3.3 Comparison with the case N = 2 

It is instructive to compare the cases > 3 and = 2 for problem (PDEm) 
with q = jj. The behaviour is the same for both in the subcritical case since 
solutions converge to a unique steady state and also in the supercritical case 
since solutions blow up in finite time. But for the critical case, the qualitative 
behaviour differs strongly. Indeed, for N — 2, blow-up occurs in infinite time 
whereas for iV > 3, there is still convergence to a regular steady state whose 
corresponding cell density has support strictly inside B, a phenomenon which 
never occurs for N = 2. We would like to suggest an "explanation" for this. 
Denoting Sm the set of stationary solutions for N > 2, we can see that we could 
as well define the critical mass as 

M = sup ||C/||oo,[o,i]- 

UGSn 

The main difference is that this supremum is not reached for N ^ 2 whereas it is 
for N > 3, which allows us in the latter case to find a supersolution that prevents 
blow-up. Thus, convergence or infinite-time blow-up seems to be determined by 
whether or not the critical mass is reached by stationary solutions. 

1.3.4 Related literature for porous medium type difTusion 

Finally, we would like to make the link between our work and the article [S] of A. 
Blanchet, J. A. Carrillo and P. Laurengot (see also the article 9 of J. Bedrossian, 
N. Rodriguez and A.L. Bcrtozzi for further results in this direction). It will allow 
to identify a formula for the critical mass M. 

The authors there study in the whole space for A^ > 3 the following Patlak- 
Keller-Segel system [PKS^) with porous-medium like nonlinear diffusion : 

= div[\JpP - ^iVc] <> a; e (22) 
-Ac = ^ t > a; e (23) 

where is the cell density and c the concentration of the chemoattractant. 
They could show that for the critical exponent p — 2 — , the system [PKS^) 
exhibits a critical mass phenomenon. More precisely, denoting m the total mass 
of the cells, they have shown the existence of Mc such that : 

• If m < A/c, solutions exist globally. 

• If TO = Mc, solutions exist globally in time. Moreover, there are infinitely 
many compactly supported stationary solutions. 

• If TO > Afc, there are solutions which blow up in finite time. 

A. Blanchet and P. Laurengot also proved in [5] the existence of self-similar 
compactly supported blowing-up solutions for to e (Mc, M+] where M+ > Mc- 
See also [38, by Y. Yao and A.L. Bertozzi for recent formal and numerical re- 
sults on self-similar and non self-similar blow-up for a generalization of system 
(PKSP) with kernel of power-law type. 

We can then observe similarities between both problem {PKS2_ ) and [PKS'^^^ ) 
for A^ > 3. Here, we would like to thank P. Laurengot for suggesting us that 
both problems should share the same stationary solutions, as we can indeed 
verify, at least formally : denoting Kq_ the Dirichlet kernel of the Laplacian in 
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a bounded domain il, it is easy to see that the steady states of (PKSq) and 
(PKS^) are respectively the solutions of 



•,1-9 



1-9 



- Kq * p = C 



and 



p-1 



where C and C are any real constants. Hence, the map p h- > /i := (2 — q) 9 p 
defines a correspondence between steady states of (PKSq) with constant C and 
steady states of {PKS'^-'i) with constant C" = (2 - g) 9 C. 
Then the formula for given in [5] also gives a formula for the critical cell 
mass M in our case : 



M = 



C,{N -I) 



or equivalently M = 



Nun' 



where C* is the optimal constant in the following variant of the Hardy-Littlewood- 
Sobolev inequality : for aU e Li(M^) n L^-f (rW)^ 



<c4h\\:_i\\h\\{ 



2/N 



Note that if we roughly put A'^ = 2 in the above formula, then we recover the 
well-known critical mass M = Stt since in this case C* = 1. 



It is interesting that, in spite of similarities in the qualitative behaviour, the two 
problems {PKSj2_) and {PRS^'-) seem to require different techniques. We 
would like to point out that our results are restricted to the radial setting, but 
on the other hand, they give a fairly more precise asymptotic description. 



Throughout the rest of the paper, we assume : 



m > 0, iV > 3 and 5 = — . 



2 Preliminaries 

In this section, we give some useful results on problem (PDEm) and two other 
auxiliary parabolic problems. For proofs, see |29] . 



2.1 Wellposedness of problem (PDEm) 

Before stating the existence and uniqueness of classical solutions for problem 
(PDEm), we need to fix some definitions. 

Definition 2.1 

Yrn = {w G C([0; 1]), u nondecreasing , u'(0) exists, m(0) = 0, — m}. 
Definition 2.2 Let T > 0. 

We say that u is a classical solution of {PDEm) on [0,T) with initial condition 
uo e Ym if •■ 

. u e C([0, T) X [0, 1]) n C\{0, T) X [0, 1]) n C^'^{{Q, T) x (0, 1]). 
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• u{0) = Uq. 

• u{t)eY,nforte [0,r). 

• u satisfies ^ on {0,T) x (0,1]. 

Definition 2.3 For any real function defined on (0, 1], we set : 

N[u] = sup — 

a:e(0,l] ^ 

Theorem 2.1 Let q e (0, 1), m > and uq E Y„i. 

i) There exists Tmax = Tmaxiuo) > and a unique maximal classical solution 
u of problem [PDEm) with initial condition uq. 
Moreover, u satisfies the following condition : 

sup Vt ||u(<)||ci([0,l]) < f°'^ °''^y ^ ^ {0,Tmax)- (24) 

te(o,T] 

ii) Blow-up alternative : Tmax = +00 or lim M[u{t)] — +00. 

Hi) Ux(t,0) > for all t G {0,Tmax)- 
Moreover, a classical comparison principle is available for problem (PDEm)- 
Lemma 2.1 Let T > 0. Assume that : 

• u^,u2 e c([o,r] X [o,i])ncH(o,r] x [0,l])nCi^2((o,T] x (0,1)). 

• For all t e (0,r], Ui{t) and U2{t) are nondecreasing. 

• There exists «o G {Ij 2} and some 7 < ^ such that 

sup f ||uio(*)llci([o.i]) < oo- 
te(o,T] 



and suppose moreover that : 

{ui)t < x'^^ir(ui)xx + ui{ui)l 

{U2)t > X^~^{u2)xx + U2{U2)1 

ui(0, x) < 1*2(0, a;) 

Ulit,0) < U2{t,0) 
Ml(i,l) < U2{t,l) 

Then ui < U2 on [0,T] x [0, 1]. 



for all {t,x) E {Q,T]x (25) 

/or an (i, a;) e (0,r] X (0,1) (26) 

for all X e [0, 1] (27) 

fort>0 (28) 

fort>0 (29) 



2.2 Two auxiliary parabolic problems 
2.2.1 Problem {tPDE„,) 

We now introduce an auxiliary transformed problem (tPDEm) which will be 
helpful in order to get some estimates implying the compactness of the trajec- 
tories in the subcritical and critical case m < M . This transformation was also 
important in [5^ in order to establish the blow-up alternative (see Theorem l2.2l 
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ii) below) , a property which will be used in the global existence part of the proof 
of Theorem O 



Denoting B the open unit ball in M^+^ and Z„i = {w E C{B), w\gB = ni}, we 
define ^° ^ ^ ^™ where | -°(^) = ^^V^^m 
We now make the following change of unknown 

^{t,y) ^ "^^i;,!^'"^ iiyeB\{0} (30) 

in problem {PDEm) with initial condition uq G Ym and obtain the following 
problem (tPDEm) with initial condition wq = 6q(uo). 

Definition 2.4 Let m > and T > 0. 

Let Uq G Y„i and wq — 9o{uq). 

We define problem (tPDEm) with initial condition wq by : 

Wt = Aw + N'^w{w + ^^)'' on{0,T]xB (31) 

w{0) = Wo (32) 

M' + ^>0 on{0,T]xB (33) 

w = m on [0, T] x dB (34) 

A classical solution on [0,T] for problem (tPDEm) with initial condition wq is 
a function 

w e C([0,r] X B)f]C^^^{{0,T] X B) 

such that (IHI|)(I211)(ISS1)(IM1) are satisfied. 

We define analogously a classical solution on [0,T). 

Let us give now the corresponding wellposedness result. 

Theorem 2.2 Let uq e Z,„ and wq = 9o{uo). 

i) There exists T* = T*{wo) > and a unique maximal classical radially 
symmetric solution w of problem (tPDE) with initial condition wq. 
Moreover, w satisfies the following condition : 

sup Vi \\w{t)\\ci{[o.i]} < oo for any T {0,T*). (35) 

te{o,T] 

ii) Blow-up alternative : T* — +oo or lim 11^(^)11^ jj — +oo. 

Connection between problems {PDE„i) and (tPDEm) ■ 
Let Wo — Ooi^o) with uo G Ym- Then, 

T^aM^ N^T*{wo). 

Moreover, for all {t,x) G [0,Trnax) x [0, 1], 

u{t,x) ^ X w{^,x^), (36) 

where we write / = /(| • |) for any radial function / on B. 
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2.2.2 Problem {PDE^^) 

We define an approximate problem (PDE^^) of problem [PDEm)- It will be 
useful since we can easily find a strict Lyapunov functional for (PDE^) and 
then prove by this way that T — lim is a strict Lyapunov functional for the 

dynamical system induced by problem [PDEm) in the subcritical and critical 

case m < M. 

Definition 2.5 Lete>0. We set : 

= (.T + forx>Q 
Observe in particular that < < for all x G [0, +oo). 

Definition 2.6 Let e > 0, m> and T > 0. 

We define problem (PDE^) with initial condition uq G Ym by : 

Ut=x^-^u,^ + uMu,) on {0,T]x [0,1] (37) 

u{0) = uo (38) 

u(<,0)=0 for allt e [0,T] (39) 

u{t, 1) = TO for all t e [0, T] (40) 

A classical solution on [0,T] of problem (PDE^) with initial condition 
Uo G Ym is a function 

u'eC{[0,T] X [0,l])f|Ci((0,r] X [0,l])f]C'^'{{0,T] X (0,1]) 
such that (inZI)(I2Hl)(IMl)(Iini) are satisfied. 

A classical solution of problem {PDE'^) on [0,r) is defined analagously. 
Theorem 2.3 Let m> 0, e > 0, K > and G with M[uo] < K. 

i) There exists T^^^ = T^^^{uo) > and a unique maximal classical solution 
u'^ on [0,T^(j^) of problem (PDE^) with initial condition uq. 
Moreover, m*^ satisfies the following condition : 

sup Vt |N/(t)||ci([o.i]) < oo for aZ/Te (0,T^,J. (41) 
*e(o,T] 

a) Blow-up alternative : T^^^^ = oo or lim Af[u'^{t)] = +oo. 

lit) {u')^>{) ont&{0,T^^^) X [0,1]. 

iv) u" e C2 iiO,T^^J X (0,1]). (not optimal) 

Connection with problem {PDE,n) ■ 
Let us fix an initial condition uq € Ym- 

The next lemma shows the convergence of maximal classical solutions u'^ of 
{PDE^) to the maximal classical solution of {PDEm) in various spaces. 
These results are essential in our proof that T = lim J^g is a strict Lyapunov 

functional in the case m < M . 
Lemma 2.2 Let uq E Y„i. 

i) Tmax(uo) < T^axM for any e>0. 



12 



a) Let [to,T] C {0,Tmax{uo))- 

a) u" — >u in C^^^{[ta, T] x (0, 1]). 

Moreover, there exists K > independent of e such that 
for all {t,x) e [to,T] x (0,1],, < 

/3) {u')^^u, mC{[tn,T] X [0,1]). 
7) {u%-^ut mC{[to,T] X [0,1]). 

3 The set of stationary solutions of problem [PDEr, 

We begin by studying the steady states of problem (PDEm), i.e. the solutions 
of the following problem : 

x^-^u" + u{uy = a; e (0, 1] (42) 
u(0) = (43) 
■u(l) = m (44) 

As is customary in evolution problems, this is essential in order to understand 
the large-time asymptotics of solutions of problem {PDEm)- 

Note that, even if we will use the results of this section only for g = they 
are all valid for any q e (0, 1). 

3.1 Existence of a steady state depending on m > 

In view of proving Theorem II. 1[ we first need to study the following Cauchy 
problem. 

Definition 3.1 For a > 0, we define the problem [Pa] by : 

a;2"«u" + u{u'f = (45) 
u(0) = (46) 
u'(0) = a (47) 

Definition 3.2 Let R>0. 

We say that u is a solution of problem (Pa) on [0] R[ if : 

• ue C\[0;R[)nC^{]0;R[). 

• u is nondecreasing. 

• u satisfies (|45l) on ]Q]R[ and also conditions (1461) and (j47D . 

This definition can obviously be adapted for the case of a closed interval [0; R\ 
or for R = +00 . 

We first state some a priori properties of the solutions of {Pa)- 

Lemma 3.1 Let a > and R > 0. 

Let u be a solution of (Pa) on [0; R[. Then : 

i) If there exists xq e [0;i?[ such that u'{xq) — then u{x) — u{xo) for all 
X G [xo; R[. 
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a) For all X G [0; R[, < u{x) < 2 and < u'{x) < a. 

Proof : i) We first note that since u' > 0, then, thanks to (|l5)) . u is concave on 
]0; R[. If xo > then for aU x e [xq; < u'{x) < u'{xo) = 0. If xq = 0, we 
have to use in addition the continuity of w' at a; = 0. 
ii) If a = then by z), u — 0, hence the resuh. 

Now, let us suppose a > and fix a; G [0; R[. u' is nonincreasing on ]0; B\ and 
continuous at x = 0, so < u'{x) < a. Hence, we just have to prove that 
u{x) < 2. 

Let us denote xq — sup{x e [0; R[, u'{x) > 0} > 0. Let x G [0; a;o[- 

First note that by z), u'{x) > 0. Since u is concave, then u'{x) < so 

u"(x) 1 ^ 

> . Hence, by integration, we get 



1 1 1 „ 

> — —— > X« 



u'{x)'i ~ u'{x)i u'(0)9 
and eventually xu' {x) < 1. Since u"{x) in nonpositive and 

uu'{x) = -xu" {x){xu' {x)Y-\ 
then uu'{x) < —xu"{x) = {u — xu')' {x). Finally, by integration, 

Uj^ ix) 

— ^ — < u{x) — xu'{x) < u{x) 
so u{x) < 2 for all x G [0; a;o[- 

If Xq = R, all is done. Else, u'{xo) = and u(xo) < 2 by continuity, so by i), 
u{x) — u{xo) < 2 for all x ^ [xq; R[. □ 

We will prove that solutions of problem (Pa) exist on [0,00). We begin by 
showing the local existence. 

Lemma 3.2 Let a > and r > 0. 

If T is small enough, there exists a unique solution of (Pa) on [0;t]. 

Proof : If a = 0, then from the previous lemma i), it is clear that is the unique 
solution of the problem (Pq) on [0, 00). 
If a > 0, let us note 

Ea = {ueC\[0,T]), u(0) =0, u'{0) = a, \\u' - a\\oo < ^}. 

Ea equipped with the metric induced by the norm ||w||_e„ = ||w'|joo is a complete 
metric space. Any u G Ea is nondecreasing on [0,t] since u' > f . 
It is clear that the following function F is well defined : 

F:Ea^C\[0;T]) 

p( V ^ r nu{s)u'{sy 

F{u)[x)=ax- / —-dsdy. 

Jo Jo s 

3 

Since for all u G i?a, ||w'||oo < -^a, we easily get that 

||W-«||oo<(^a)^+^^<| 
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if T is chosen small enough. Hence, F sends Ea into Ea- 
By the mean value theorem, since for all u ^ Ea, ^ < u' < 



Hence, if r is small enough, is a contraction so there exists a fixed point of 
F. Since F{u) G C^{{0,t]) when u e C^{[0;t]) then m is a solution of (Pa). 
Eventually, it is easy to check that a solution is necessarily a fixed point of F, 
which proves the uniqueness. □ 

Remark 3.1 Let u the local solution of {Pa) on [0,t]. We would like to stress 
that u is not up to x = 0. Indeed, one can see that 

u'{x) =a x'^ + oix'^). 

q 

Moreover, one can obtain an expansion of u at any order in powers of . We 
proved in f29\l that solutions of problem (PDEm) share these properties with the 
stationary solutions. 

This fact is important because, since the bacteria density p is proportional to 
Ux, it explains the "emaciated" shape of p that can be observed on numerical 
simulations. 

Theorem 3.1 Leta>0. 

There exists a unique maximal solution of (Pa). 
Moreover, it is globally defined on [0,oo). 

Proof : The case a = is obvious. 

Let a > 0. By the previous lemma, for a given r small enough, we have a 
unique classical solution Ur of {Pa) on [0; r]. We can now consider the following 

T 

ordinary differential equation on the interval [— , +cx)[ : 

u 



x^-1 

u{^) = M^) 
u'{^)^K{^)>o 

Setting W — {u, u'), it becomes a first order nonautonomous ordinary differen- 
tial equation 

W' = F{x, W) 

where F{x, u, v) = {v, — - — ). 

Let us denote il — Rx]0;+cx)[. Since F is locally Lipschitz continuous with 
respect to W in il, the classical Cauchy-Lipschitz theory can be applied. Then, 
there exists a maximal solution 

ueC'{[^;X*)) 
of the ordinary differential equation such that 

u{l)^Ur{l), U'C-)=U'^{^) 
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and 

{u{x),u'{x)) e n for all x e [^;^*). 

By local uniqueness in the classical Cauchy-Lipschitz theory, Ur — u around 
a; = so we can glue Ut and u and get a solution of problem (Pa) on [0; X*). 
If X* = +00, all is done. So we suppose that X* < 00. 
Since u in nondecreasing and bounded above by Lemma l3. II ii). then 

/ = lim u(x) exists. 

Hence, we can extend continuously the function u by setting u{x) = I for x > 
X*. But {u{x),u'{x)) must leave any compact of 51 as a; goes to X* so, by 
Lemma [5TT|) ii), the only possibility is that lim u'{x) — 0. So u e C^{[0', +00)). 

And now, thanks to lim u"{x) — 0, so u £ C'^((0;+oo). Moreover, u 

£C— fX* 

clearly satisfies (|45p on (0, +00) and is then a global solution of (Pa)- 
Let V another global solution of (Pa)- By the result of uniqueness around x = 
and the uniqueness due to classical Cauchy-Lipschitz theory in 57, u and v 
coincide for all a; e [0,X*). li X* —00, all is done. Now, assume that X* < 
+(X). v{X*) = u{X*) by continuity. As v'{X*) = u'{X*) ^ 0, then by Lemma 
EH) i), v{x) = v{X*) = u{X*) for all x > X*. Hence, v = u. □ 

Notation 3.1 Leta>0. 

We denote Ua the unique solution of (Pa) on [0,+c»). 
Lemma 3.3 Let a > 0. 

There exists xq >0 such that Ua{x) = Uaixo) for all x > xq. 

Proof : Suppose the contrary. By Lemma it implies that U^{x) > for 
all a; > 0. Since Ua is nondecreasing and has 2 as an upper bound, there exists 
I < 2 such that Ua{x) tends to / as a: goes to infinity. 

As U'a is nonnegative and nonincreasing, then J7^(a;) has a nonnegative limit as 
X goes to +00, but this limit has to be since Ua is bounded from above. 
Moreover, for all a; > 0, 

4-U'a{-f-''^Ua{x)-'' ' 



dx °' ' ' 'dxx^-1' 

Let Xq > and a; > a;o. By monotonicity, Ua{xo) < Ua{x) < Z so by integration 
on [a:o;a;], 

^ - ^) < Ki-f-' - Ki-oY-' < - ^) ■ 

Finally, we let x go to ~\-oo and then go to +00 to obtain that 

1 

I— 



KM 



a;o->+oo a;o 

Then, Ua{x) goes to infinity as x goes to oo, which is a contradiction. □ 

Proposition 3.1 Ua{x) — Ui{ax) for all x > and all a>0. 

Proof: Let V{x) = Ui{ax). Clearly, V G Ci([0, +oo))nC2((0; +00)), ^(0) = 0, 
y (0) = a and for a; > 0, 

x^-^Vix) = {axf-m['{ax)a'i = f/i(aa;)J7{ (aa;)«a« = Vix)V'{xy. 

The result then follows from the uniqueness of the solution of problem (Pa)- O 
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Remark 3.2 Behind this proof is the fact that La and D commute, where for 
any u £ ([0, cxd)), a > and x > 0, La{u){x) ~ u{ax) and Du = xu^- 

This proposition drives us to the natural fohowing definition. 

Definition 3.3 We call M the critical mass. 

Note that M is the maximal value of each Ua, for all a > 0. 

We also will use the following notation. 

Notation 3.2 We denote by A the first x > such that Ui{x) — M. 
Proposition 3.2 

i) If m Cz [0;M), there exists a unique a € [0; A) such that C/a(l) = w. 

'ii) Ua{l) = M if and only if a > A. 

Proof : i) C/q(1) = Ui{a) and Ui is a bijection from [0; A) to [0; Af). 

ii) t/a(l) = M if and only if Ui{a) = Af, which is equivalent to a > A. □ 

From this follows Theorem 11.11 announced in the introduction. 

3.2 Order and topological properties of the set of station- 
ary solutions 

Now, we shall describe two simple but very important properties of the set of 
stationary solutions : it has a total order and its topology behaves well with 
respect to that order. 

Notation 3.3 Letm>0. 

• — {u £ ^"^([0; 1]), u nondecreasing , u{0) — 0, u(l) = m} is the com- 
plete metric space equipped with the distance induced by the norm. 

• Vm{u, e) C Yrn IS the open ball of center u G Ym and radius e > 0. 

• We say that two functions m, u G Ym satisfy u < v if u < v and u ^ v. 
If u £ Yra and V C Ym, we say that u <V if for all v Cz V , u < v. 
Similarly, we define V < u. 

Proposition 3.3 

i) Suppose < a <b. Then Ua < Ut. 

The set of stationary solutions is a totally ordered set. 

ii) Suppose A < a < b. 

a) There exists e > such that VM{Ua,e) < Ut. 

j3) There exists e > such that {u G YIj, u < Ua} H VniUb, e) = 0. 

Proof : i) Let x>0. Ua{x) = Ui{ax) < Ui{bx) = Ub{x), so Ua < Ut. 
Since, U^{0) < J7^(0), then Ua < Ub- 

ii) a) There exists 7 > such that for aU x G [0; 7], ^^x < Ub{x). 
Let us set ei = £2 = ^ niin [Ub{x) — Ua{x)] and e ~ min(ei,e2)- 

Note that e > since for all x G [7; ], 

Ua{x) = Ui{ax) < Ui{bx) = Ub{x) 
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because Ui is increasing on [O;^]. 

Let u e VM{Ua,e)- Since u e Va,/ (f/a, ei), then 

, , b ~ a b — a a + b 

Hence, for all x S [0;7], 

, , a + b rr / \ 

u(x) < ^ X < Ub{x). 



Since u € VM{Ua, £2), it is clear that for all x £ [7; u(x) < Ub{x). Moreover, 

lence, 

u < Ub- 



b - 

if x > f then u{x) < M = Ub{x). Hence, 



13) Let u GYm such that u < Ua- Then w'(0) < a. Let us set e = ^ > 0. 
Since \\Ub - w||ci > Ul{Q) - u'(0) >b-a then u i VM{Ub,e)- □ 



4 Convergence to a stationary state in critical 
and subcritical case m < M 

All this section only concerns the case m < M . We recall : 

Definition 4.1 — P| C^([0, 1]). Equipped with the distance induced by 
the norm, Y^ is a complete metric space. 

We shall prove that problem {PDE,n) defines a continuous dynamical system 
on Y^ which admits a strict Lyapunov functional. We shall be able to prove 
that classical solutions of (PDEm) converge to a stationary state as times goes 
to infinity, even in the case m = M where there is a continuum of steady states. 



4.1 Estimates 

Lemma 4.1 Let m < M and uq E Y„i. Then, 

Tmax{uo) = 00 • 

Moreover, for each K > 0, there exists Ck > such that for any uq G Ym with 
Af[uo] < K then, 

sup N[u{t)] < Ck- 

tG[0,oo) 

Proof : Let T,nax = Tmaxiuo). 

From Theorem 12. 1[ it is sufficient to prove that 

sup J\f[u{t)] < 00. 

te[o,T„„,) 

This fact easily follows from a comparison with a supersolution of problem 

{PDEm)- The main idea is that since m < M, if a is large enough then uq < Ua 

and Ua is then a supersolution so < u{t) < Ua for all t e [0, T,nax)- 

More precisely, since uq is differentiable at a; = 0, x t-^ i^oX^) extended 

continuously to [0; 1], so there exists K > M such that A/'[uo] < K. 

Let us set a = -j^^. 
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For X £ [0; uo{x) < Kx < Ua{x) since by concavity, Ua is above its chord 
between x = and a; = ^. On 1], uq < m < M ^ Ua- 
Hence, uq < Ua on [0, 1]. Finally, since M > m, C/q is a supersolution of the 
PDE, so u{t) < Ua for all t S [0,Tmax)- By concavity of Ua, we see that 
J^[Ua] = a, so sup Af[u{t)] < a < oo. Then T„iax ~ oo. 

We notice that the choice of a depends only on K, whence the second part of 
the lemma. □ 

Before going further, we would like to recall some notation and properties of 
the heat semigroup. For reference, see for instance the book of A. Lunardi. 

Notation 4.1 

• B denotes the open unit ball in R^+^ . 

• Xo = {We C(B), W\aB = 0}. 

• {S(t))t>o denotes the heat semigroup on Xq. 

• (X0)eg[Q.i] denotes the scale of interpolation spaces for {S{t))t>o- 
Properties 4.1 

. Xi ={WeC\B), W\dB=0}. 

• Let 70 e (0; i]. For any 7 e [0, 270), 

with continuous embedding. 

• There exists Co > 1 such that for any 9 £ [0; 1], W £ C{B) and t > Q, 

\\S{t)W\\xe<^\\WU. 

We just want to introduce a specific notation we are going to use. 

Notation 4.2 Let {a,b) G (0, 1)^ We denote L{a,b) = j^^^- 
For all t > 0, /„* j^^^ = t^—H{a, b). 

We will now give an estimate from which follows a compactness result. 

Lemma 4.2 Let m < M, 7 e [0; 1), to > and K > 0. 

There exists Dx > such that for any uq G Ym with Af[uo\ < K , then 

sup|iw(t)|ipi,^ <Dk. 

t>to 

As a consequence, {u{t), Af[uo] < K, t> to] is relatively compact in Y^. 

Proof : Let uq G Y^ such that A/'[uo] < K. Let wq = 9q{uq). 
First step : from Lemma 14.11 there exists Ck > such that 

sup Af[u{t)] < Ck. 

tG[0,oo) 
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Since for t>0, \\w{t)\\^g ~ JV[u{^)], we deduce that w is global and that 

sup sup J^[u{t)] < Ck- 

te[o,oo) ' te[o,oo) 

Second step : Let t — and t e [0, r]. 
Denoting Wq = wq — m, then 

w{t) -m^ S{t)Wa+ S{t~s)N^w[w + ——] ds, (48) 



so 



TV 



<m+^{CK + m) + N^ I -^^CK\\w{s)\\''^,ds. 
V r Jo Vt ~ s 

Setting h{t) — sup Y^||w;(s)||pi, we have h{t) < oo by (1551) ^-^id 
se(o,t] 

/■* 1 

Jo s-i^/t — s 

^^t\\w{t)\\c^ < mV^+Coim + Ck) + N^CkCdI Q, |^ t^-ih{ty. 
Let T e (0,r]. Then, 

h{T) < + CD(m + Ck) + iV^C^Cz^/ Q, |^ T^-'-hiTf. (49) 

Setting yl = m^r + Cc(m + C^) and B = N^CkCdI (5, f ) 2^, assume that 
there exists T G [0, r] such that h(T) = 2A. Then, 



^1-9 < 5^1-1 which imphes T > 



B 



Let us set r' = min f r, | ^ 



Since /i > is nondecreasing, /lo = lim h{t) exists and ho < Ahy (|49l) . So by 

t-»-o+ 

continuity of h on (0,t'], /i(t) < 2 A for ah t e (0,t'], that is to say : 

2A 

\\w{t)\\ci < —j= for ah t e (0,t'], 
\Jt 

where A and r' only depend on K . Then, setting Ak — 2A, we have 

sup Vi\\w(t)\\ci < Ak. 
te[o,T'] 

Third step : Let 70 G (^, i) and t e [0,t']. 

Setting W = w — m and Wq = — to, then for f > 0, due to (051), we get 



^■D , ^ , Ar2 /" Cd ^ {Ak 



Then we deduce that : 

ti+^°\\W{t)\\x. <CDiCK + m)+N^CKCDiAK)Hi+-"' f ^- -ds 

< Coim + Ck) + N^CKCoiAKyH^ + 70, |)t'^'-^\ 
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Hence, since ^^+-y„ C C^ ''{B), we deduce that there exists A'^ > depending 

only on K such that sup t^^'' \\w{t)\\ fji^^ < A'j^. Then, 

te[Q,T'] 



A' 



Last step : Let t' > j^. Since t' < we can apply the same arguments by 
taking wq it' — r) as initial data instead of wq , so we obtain 

for all t'>^, imOlici,.(U) 

Eventually, coming back to u(t), thanks to formula p6p . we get an upper bound 
Dk for ||M(t)||pi,-x valid for any uq e Ym such that A/'[uo] < K. □ 

4.2 A continuous dynamical system {T{t))t>o 

We recall the definition of a continuous dynamical system on Y^. 
For reference, see [151 chap. 9, p. 142]. 

Definition 4.2 A continuous dynamical system on Y^ is a one-parameter fam- 
ily of mappings {T{t))t>a from Y^l^ to Y^ such that : 

i) T(0) = Id. 

it) Tit + s)= T{t)T{s) for any t,s>0. 
Ill) For any t > 0, T{t) G C(r„\, F^). 
iv) For any uq G Y„,, 1 T{t)uo G C((0, cx)), F^). 

Remark 4.1 Continuity at t — is sometimes included in the definition, but 
it is not required for our needs. 

Definition 4.3 Let uq G Y^ and t > 0. 

We define T{t)uo — u{t) where we recall that u is the classical solution of 
problem (PDEm) with initial condition uq. 

Proposition 4.1 {T{t))t>o is a continuous dynamical system on Y^^. 

Proof : Thanks to Lemma [4.11 we know that T{t) is well defined for alH > 
and by definition of a classical solution, Tit) maps Y^ into Y^. 
ii) is clear by uniqueness of the global classical solution, 
iv) comes from the fact that u G C((0, oo), C^([0, 1])). 
in) Let t > 0, uo G and (u„)„>i G Y^. 

Assume that m„ — >■ uq. Let us show that u„(i) — >■ u{t). 

n— foo n— >-oo 

We proceed in two steps. 

First step : We show that if u„ — > uq, then u„(i) — > u{t). Let rj > 0. 

By (IMl), there exists C > such that for all s G (0,t], ||Ma;(s)||oo < 
So we can choose 77' > such that 

^'e/o[il'"-(^)llL+i] ds ^ ^_ 
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Let no > 1 such that for aU n > uq, ||u„ — uo||oo,[o,i] ^ 

Let n > riQ and s G [0,t]. We denote Un{s) the solution at tune s of problem 
{PDEm) with initial condition m„ and set : 

z{s) = [u„(s) - w(s)]e-/o[ll''-(''')ll~+il . 

We see that z satisfies 

Zs = x^^~ z^x + bzx + cz (50) 

where 

b{s,x) = Un{s) (""-^^('^'^-^ "'^^'^'^1 if (u„)^(s,a;) ^ U2;(s,a;) and else 

(Un)x{S,x) - Uoc[S,X) 

and 

c=[K)^-|lzi.||^-l]<0. 

Since z G C([0, i] x [0, 1]), z reaches its maximum and its minimum. 
Assume that this maximum is greater than Since z = for x = and x = I 
and z < 77' for s = 0, it can be reached only in (0, t] x (0, 1) but this is impossible 
because c < and (jSOp . We make a similar reasoning for the minimum. Hence, 
\z\ < rj' on [0,t] X [0,1]. 

Eventually, - ?/|U,[o,i]x[o,*] < r/'e/o < rj for aU n > tiq. 

Whence the result. 

Second step : since u„ — > uq, ||u„||(71 is bounded so there exists K > 

n— f 00 

such that for all n > 1, Af[un] < K. Then, from Lemma 14.21 since t > 0, 

{un{t), n > 1} is relatively compact in and has a single accumulation point 

u{t) from first step. Whence the result. □ 

4.3 A strict Lyapunov functional for {T(t))t>o 
4.3.1 Reminder on strict Lyapunov functionals 

We recall some definitions in the context of a continuous dynamical system 
{T{t))t>o on y^, including strict Lyapunov functional and Lasalle's invariance 
principle. 

Definition 4.4 Let uq G F^. 

• 71 (wq) = {T(t)ua, t > 1} is the trajectory of uq from t = 1. 

• uj{uq) ^ {v e F„\, 3tn -> +00, t„ > 1, T{t„)uo — > V in Y^} 

n— f +00 

is the uj-limit set ofuQ. 
Definition 4.5 

i) J- (z C{Y^,M.) is a Lyapunov functional if for all uq G y^, 
1 1— )• T[T{t)uo\ is nonincreasing on [1,+cx)). 

ii) A Lyapunov functional T is a strict Lyapunov functional if 

J-[r{t)uQ\ — J^[mo] for all t > implies that uq is an equilibrium point. 



22 



Proposition 4.2 Lasalle's invariance principle. 

Let Ufj G Y^. Assume that the dynamical system (T{t))t>o admits a strict 
Lyapunov functional and that 71 (mq) is relatively compact in y„\. 
Then the co-limit set w(uo) is nonempty and consists of equilibria of the dynam- 
ical system. 

See [m p. 143] for a proof. 

4.3.2 Approximate Lyapunov functionals 

We recall that for all e > and a; > 0, fe{x) — {x + e)'^ — e"^ . In order to introduce 
the approximate Lyapunov functional, we first need a double primitive of 

2,2-1? 

■J- that converges uniformly to H{x) = -7 ^ on compacts of M+. The 

3' _ ^ {2-q){l~q) 

next lemma provides it. 

Lemma 4.3 Let q e (0, 1) and e G (0, 1]. 

For x>0, we set H^{x) = -f7{t)^y + 1^ ""'^ ^(^) = (2 ~ q){l - q) ' 

i) is continuous on [0;+oo), twice differentiable on (0;+oo). 
H'^ = 77 on (0,+oo). 

ii) converges uniformly to H on [0; i?] as e tends to 0, for any R > 0. 

Proof : Let i? > and x E [0, R]. We begin with two remarks : 

- Since f^ > fi and /i is concave, then, denoting K — we have 

ft{t)>Kt foranyt>0. 

- We also note that for any i > 0, 

0<i«-/.(i)<e'. 
Indeed, setting g{e) — {t + e.y — e"^, we have 

- = g(0) - .9(6) = g\es)ds = e C - —^-—ds < e'^ . 



(se)i-« (t + es)i- 

i) Let us set for y > 0, ^e{y) = Ji j^- Since for all t > 0, f^{t) > Kt, then 
7(:(y) = 0(|log(y)|). Hence 7^ is integrable on {0,R]. Then, is continuous 
on [0, i?]. The other facts are obvious. 

ii) We can write 

Using the first two remarks, we get 

Whence the result. □ 
Definition 4.6 Let e > 0. For u G Y^, we define 



Te{u)=l H,{u:,) - ^-^dx. 
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Lemma 4.4 Let uq eY^. For all < t < s < T^^^ 
More precisely, for all t > Q, 



dx. 



Proof : Let t > and 77 > such that / = [t - ry, i + 77] C (0, T^,ax)- 

By Theorem 12.31 iii). there exists /i > such that {u'^)x > /i on / x [0; 1] so it 

will allow us to use that £ C^([/i, +cx))). 

By Lemma [2.21 {u'^)tx, (w'^)t, {u'^)x are bounded on / x [0; 1] and there exists 
K>0 such that \iu')^,\ < ^ on / X [0; 1]. Moreover, < ^ < ^^^M^ < 
on / X [0; 1] for some K' > 0. Note also that ul^ — u^^ since G 
C^{{0,T) X (0, 1]). All these facts allow us to differentiate and then to 

integrate by parts : 



t 



d 

Hence, -T.iu'it)) J 



^ since {u%{t, 0) = {u%[t, 1) = 0. □ 



x^-if,{{u%) 



4.3.3 A strict Lyapunov functional 
Definition 4.7 For u G Y^, we define 



„2-9 

where H(v) = for all u € M. 

(2-q)(l-g) 

Theorem 4.1 

T is a strict Lyapunov functional for the dynamical system T{t)t>o on Y^. 

This fact cannot be obtained directly by the formal computation shown in the 
introduction since can vanish on a whole interval. Indeed, consider for in- 
stance an initial condition uq e Ym such that uq > Ua where a > A. By 
comparison principle, for all t > Q, M > u{t) > Ua so Ux{t) = at least on 

La' J 

Proof : First, by application of the dominated convergence theorem, since H is 
continuous on [0,oo), it is clear that is continuous on Y^^. 
Secondly, let Uq £ Y^ and < i < s. Let us show that F{u{t)) > J^{u{s)). 
Let t' > 0. We first need to show that — > T{u{t')). 

Lemma[221tellsusthat'u'(t') — > u{t') inCi([0,l]). Li particular, (u<^(i'))ee(o i) 

e— >0 ' 

is bounded in C^([0; 1]) so we have a domination independent of e. Since 
converges uniformly to H on compact subsets of [0; +00) and H is continu- 
ous on [0; -|-oo), by the dominated convergence theorem, we obtain easily that 
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e— fO 

Then, as we know from Lemma that J'g[u'^{t)] > J^g[u'^{s)], the resuh fohows 
by letting e go to zero. 

Thirdly, denoting R — sup ||7io(t')||(7i; we want to show that 

t'e[t;s] 

J^[u{t)] - :F[uo{s)] > ^ 1 1 {uo)l (51) 

[H + 1)1 J J[t:s]x[0;l] 

By Lemma [4.41 

J Jlt;s]xlO-l] ''h{{u')x) 

By Lemma {u'^)x tends to Ux in C{[t, s] x [0; 1]), so there exists eo > such 
that for all e G (0;eo), sup \\u'{t')\\ci < R + l- 

Note that ^ > 1 for < x < 1 and that < < < (R + 1)« on 

[t, s] X [0, if. So, 

[R + 1)9 J J[t;s]x[0;l] 

By Lemma [2?2l {u'^)t tends to (wo)t in C{[t; s] x [0, 1]), hence by taking the limit 

as e goes to 0, we obtain the result. 

Eventually, assume that J^[u{t)] = J^[uo] for all t > 0. 

Let [t,T] C (0,oo). Formula shows that ut = on [i,r]. So = on 
(0, oo) X [0, 1] . Then, by continuity of u on [0, oo) x [0, 1] , we get that T{t)uQ = uq 
for all t > 0, i.e. uq is an equilibrium of the dynamical system {T{t))t>o- O 

4.4 Convergence to a stationary state for m < M : proof 
of Theorem 11.21 

In the case of m < M, there is a unique stationary solution for problem {PDE„i) 
so the convergence is not really surprising. 

But for m = M there is a continuum of stationary solutions (all C^a |[o,i] for 
a > A) and the behaviour could be much more complicated. However, thanks 
to the good properties of the set of steady states (see Proposition [5T31) and since 
the problem is one-dimensional, convergence can be shown by arguments in the 
spirit of p7] or [39]. 

Proof of Theorem \1.2\ : Let uo £ Ym- Let us set ui — uo{l) e Y^. 
To get the result, we just have to study lim T{t)ui. 

t— f-f oo 

Thanks to Lemma l4.2l 71 (ui) is relatively compact in Y^ and since is a strict 
Lyapunov functional for (T(i))t>o, we know by Lemma 14.21 that the cj-limit set 
ui{ui) is non empty and contains only stationary solutions. 
First case : m < M. Then from Theorem ll.il there exists a unique stationary 
solution Ua with a < M. Hence, w(ui) = {Ua} so T{t)ui — > Ua ■ 

Second case : m — M. Assume by contradiction that io{ui) contains two dif- 
ferent stationary solutions Ua and Ub with A < a < b. Then we chose any 
c € (a, b). From Proposition l3.3p ii), there exists e > such that Va/(C/q, e) < Uc 
and {u e Ym, u<Uc}n VM{Ub, e) = 0. 

Since Ua G w(wi), there exists ta such that T{ta)ui G VMiUa,e)- Hence, 
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T(ta)ui < Uc and then by comparison principle, for all t > ta, T[t)ui < Uc- 
But, since Ut G uj{ui), there exists tb > ta such that T{tt)ui e Vjvf(C^fc,e), and 
this is a contradiction because {u G Ym, u < Uc}r\VM{Ub, e) ~ 0. Hence, uj{ui) 
is a singleton {Ua} with a > ^ so T{t)ui — > Ua- □ 

5 Finite time blow-up and self-similar solutions 
in supercritical case m > M 

In this section, we only consider the supercritical case, i.e. when m > M. 
We shall prove that classical solutions of problem (PDEm) blow up in finite 
time. The idea of the proof is to exhibit a subsolution u{t, x) — V{a{t)x) which 
turns out to be a self-similar solution after some time. 

This is why we are interested in the following ordinary differential equation. 

5.1 An auxiliary ordinary differential equation 

Definition 5.1 Let e > 0. 

We define the problem (Q^) by : 

x^-^V + VV"^ = exV x>Q (52) 
V{Q) = (53) 
V{Q) = 1 (54) 

Definition 5.2 Let e > and R> 0. 

We say that V is a solution of problem {Qe) on [0; R[ if : 

• Ci([0;i?[)nC2(]0;i?[). 

• V is nondecreasing. 

• V satisfies (j52p on ]0;i?[ and the conditions (1531) and (j54[) . 

This definition can obviously be adapted for the case of a closed interval [0; R] 
or for R ~ +oo. 

We summarize in the following theorem some very helpful results about solutions 
of problem (Qe)- 

Recall that C/i is the solution of problem (Pi) and that A is the first point from 
which Ui is constant. 

Theorem 5.1 There exists a unique solution of problem {Qe) on [0, oo). 

// e > is small enough, is concave and there is a first point < oo from 

which Ve is constant with value greater than M . 

Moreover : 

i) \\Ve - C/i||ci([o,oo)) ^ 0. 
a) Ae — > A. 

e-i-O 

Hi) VJA + 1) — > M. 

That is to say that the constant reached by can be chosen as close to 
M as we wish provided that e is small enough. 
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The proof of Theorem 15.11 follows from the successive lemmas we are going to 
state now. Note that the proof that > M and that is concave for small 
e > is included in the proof of Theorem 11.41 

We begin by giving an a priori property of the solutions of problem (Qc). 

Lemma 5.1 Let e > and V a solution of problem (Qe) on [0; +oo[. 
If for some xq > 0, V{xo) = 0, then for all x > xq, V(x) = V{xq). 

Proof : First, a;o > since V{Q) = 1, then V{x{)) > since V is nondecreasing. 
Let R — sup{r > xq, s.t. V constant on [a:o,r]}. 

Assume by contradiction that R < oo. Then, by continuity, V{R) — V{xo) > 
and V{R) — 0. Writing equation ((5^ as 

x^-iy ^ V'^i-V + exV^~''), 

we see that the first factor of the RHS is nonnegative and that the second one 
keeps negative for x close enough to R by continuity since —V{R)+€RV{RY~'^ < 
0. Hence V is nonpositive for x close enough to R. But V{R) — and V > 0, 
then V = near of R, which contradicts the definition of R. So, R — oo. □ 

We now prove the local existence of a solution of problem (Qe). 
Lemma 5.2 Let e E (0, 1]. 

There exists d > independent of e such that the problem {Qe) admits a unique 
solution on [0, S\ . 

Proof : The method used is a fixed point argument, as for the local existence 
of the solutions of problem [Pa]- 
Let us define 

E^{V eCWQ,5]), t/>0, F(0)=0, t/(0) = l, ||l>-l||oo < ^}. 

E equipped with the metric induced by the norm HI^Hb ~ Halloo is a complete 
metric space. We define F by : 

F ■.E^C^{[{y,5]) 

F{V){x) =x- r r ^^^ds dy + e r r -^ds dy. 

Jo Jo S 1 Jo Jo fil 1 

Since for ah V e E, \\V\\oc < | and < e < 1, we easily get that 

iiF(v)'-iiu<((|r' + |):^<i, 

provided that S is chosen small enough. Hence, F sends E into E. 

We can apply the mean value theorem to function z z*, since for all V £ E, 

i < 1/ < |. Eventually, we obtain for all (Vl, V^2) e E'^ : 

\\FiV2) - F{V^)\\E < f ^ + -t|— + -) S^V2 - V,\\e. 

Hence, if 6 is small enough, is a contraction so there exists a fixed point V of 
F. Since F{V) G C^{{0, 6]) when V G Ci([0; r]) then y is a solution of (Q,). 
Eventually, it is easy to check that a solution of (Qe) is necessarily a fixed point 
of F, which proves the uniqueness. □ 
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Lemma 5.3 Let e G (0, 1]. 

There exists a unique > S and a unique maximal solution Ve on [0, A^) of 
(Qe) such that > on [0,^^). 

Proof : By the previous lemma, we have a unique classical solution Vq of (Qe) 
on [0;5]. We can now consider the following ordinary differential equation on 

the interval • 

VV'i e ■ 
V = -^ + -^V 

Setting Z = (V, V), it becomes a first order nonautonomous ordinary differential 
equation 

Z = F{x,Z, e) 

where F(x, V, W,e) = { W, + — W 

\ x^~i x^~i 

Let us denote O = M x (0; +oo). Since F is locally Lipschitz with respect to Z 
in f2, the classical theory of Cauchy-Lipschitz can be applied. Then there exists 

a maximal solution of the ordinary differential equation V G C^([-,Ae)) such 

that V{^) = Vo{^), V{^) = Vo{^) and {V{x),V{x)) e Q for all x G [^,A,). 

The local uniqueness given by the classical Cauchy-Lipschitz theory ensures that 
we can connect Vq and V without any problem and then get a solution of our 
ordinary differential equation on [0; ^e). □ 

Lemma 5.4 Let e G (0, 1]. 

i) We have the following formula : for x £ [0, ^J; 

V^{x) = exp(^ — I 1 - (1 - 9) / exp( — ^ — r^)^« 



a) For all X G [0,A^), 



< V,{x) < exp(^^^-^), 



0< T' 
Jo 



0<V,{x)< I exp(^L«)£!)ds. 



Proof : i) Let us set w = V^. On (0, A^), w satisfies 

w = -z — H — z — w. 

g2-q 3.1-9 

Wc recognize a Bernoulli type ordinary differential equation. Since u> > on 
[0, Ae), then we can divide by w'^. Setting z = we obtain 

e(l-g) Ve{s){l-q) 
z — z — 

j.1-9 „2-9 
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which can be easily integrated. Whence the formula. 

ii) Ve is increasing and Ve(0) = so Ve > on [0, A^) whence the results. □ 

Remark 5.1 A consequence of last lemma is that 
tfl-{l-q) /; ^ exp(-^ii^)ds < 0, then A, < x. 

Lemma 5.5 Let e G (0, 1]. 

There exists a unique solution Vc of problem (Qe) defined on [0,+oo). 
Proof : li Ag — oo then all is already done. 

Else, if < oo then Ve{x) has to go to zero as x goes to A^. Indeed, 
{Ve{x),Ve{x)) must go out of any compact of K x (0, +cx)) as x goes to A^ 
but by ii) of the above lemma Ve{x) and V^{x) keep bounded for x bounded. 
SoV;(a;) 0. 

x—yAc 

Hence, by Cauchy criterion, Ve{x) has a limit as x goes to A^. Moreover, by 
the equation (|5^ . Ve{x) goes to zero as x goes to A^. Since the constants are 
solutions of ([5^ . then can be extended by the constant on [A^, +oo) to a 

function on (0, +00) which is a solution of problem ([5^ . 
For proving the uniqueness, let V a solution on [0, +00). By uniqueness of the 
solution on [0,A^) then V = on [0, A^). Then, by continuity, V{Ae) ^ V^iA^) 
and V{A^) = 0. But now by Lemma TS. 11 V is constant on [Ag,oo) so V = 
on [0,00). □ 

Lemma 5.6 

i) A^ < 00 for e small enough. 
Moreover A^ — >■ A. 

ii) ||Ve - J7i||ci([o,oo)) ^ 0. 

As a consequence, V^{A + 1) — > M . 

Proof : 

First step : We show that there exists 7 e (0,(5] independent of e such that 

\\y^-Ui\\cH[o..,\)^,^- 

Let us set 7 = min((5, where 5 is a short existence time for all and 
(5' = C/i > on [0, A) so equation (pS)) can be written on [0,A) as 

d Ui^'' Ui 



dx 1 — q x'^ 1 
Hence, we have the following formula 

Ui{x)^(^l-{l^q)£^dsy~\ (55) 

which is valid for all a; € [0, A], by continuity. 
Let X e [0, 7]. 

V,ix)-Ui{x) = 

(».-f--i)(i-a-„rigle.p,-i<i-i)i!M, ■ 
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+ l-(l-g) 



l-(l-9) 



Then, 



-(1-9)^" 



\V,{x) - Uiix)\ < (e^— 



1) 



=2-9 



because \a^-i - 6i-<j | < ]4^|a - 5| for all (a, 6) e [0, 1]^ since > 1 

|K(s)-C/i(s)| 



, (i~a)VL 



|V;(x)-C/i(a;)| < (e^^^ -1) 
'■^f/i(s) 



=2-9 



(1 — e 9 )as. 



Then, we get 



liV;-C/i||oo,[o,,] < (e' 
+ llV;-C/i||oo,[o,7]— + (i-e" 
Since ^ < i and < ^, then 

liv;-c>illoo,[o,7]<2 



, (i-<;)7'' 



1) 



,s2-9 



ds. 



(l-g)^'' 
e 9 — 1 



Hence, ||Ve 
Second step : Let A' < A. 

Let us show that for e small enough A^ > A' and that 

llK-C/i||cH[o,A'])^0. 

Let us denote Va,b the solution of (|5^ such that 14,6 (7) = « and 14,6(7) = ^• 
Ui is the solution of equation ([5^ for e = and initial condition (a, b) — 
{Ui{j),Ui{'-f)). Since C/i > on [7, A'], the classical Cauchy-Lipschitz theory 
is here available, and by continuity of the solutions on [7, A'] with respect to 
the parameter e and the initial condition (a, &), we know that if e is small 
enough and if the initial condition (a, &) is close enough to {Ui{'-f),Ui{j)) then 
1114. 6 — t^i||ci([7,A']) is as small as we wish. But, thanks to the first step, taking e 
even smaller, {Vt{'-f), Ve{-f)) can be made as close as necessary of {Ui{j), Ui{'-f)). 
Eventually, \\Ve — J/i||ci([o,A']) is as small as necessary when e is small enough. 
This implies of course that > A' for e small enough. Whence the results. 
Third step : Let us show that ||Ve — C^i ||ci([o,yi+i]) — ^ 0. 
Let a > and t] = (f )^"''. 
By formula ((SSl) . we have {1 — q) ds = 

there exists A! < A such that 1 — (1 — 9) ^^^ds < 77 and \\Ui 
Let X e [A', min(Ae, A + 1)]. We know that 



ds = 1 and since J7i = on [A, 00), 

00, [A' . 



V^{x) = cxp( 



£(l-g)£« - 



)(l-(l-g)/;^exp( 



)ds 



1 

1-9 
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then !> (x) < exp(^ii^^M±l)l) (l - (1 - q) ^ds) since V, > 0. 
There exists ei > such that if < e < ei, exp( ) < 2, so Ve{x) < 
2rjT^ = |. Then || V;||oo,[A',a+i] < f since V; = in [A„oo]. Eventually, if 
< e < ei, then \\Vf — t^i||oo,[yi'.yi+i] + § = From the second step, let 
£2 > such that HV'c — [/i||oo.[o.A'] < a for < e < £2- Let us set eo = min(ei, £2)- 
Then, if if < e < eo, then \\V^ - C/i|U,[o,A+i] < a. 

Fourth step : Let us show that for e small enough, < 00 and even that 
— > A. 

e->-0 

Let A' E {A, ^ + 1]. Since, by formula (|55p . we have 



and since 



l!^.-f/i||ci([o,A']) ^0, 



then there exists eg > such that for < e < eo, 

Jo s 1 

Let < e < eo- From remark [5TT1 then A^ < A' < 00. Combined with step 2, 
this proves that — > A 

Last step : Let e > small enough so that A^ < A + 1. Then, Ve — Ui = on 
[A + 1,00) thanks to Lemma f5.ll Moreover, thanks to step 3, Ve{A + 1) — > 

UiiA + 1) = M. Whence \\V, - i7i||ci([o,oo)) ^ 0. '^D 

5.2 Blow-up and existence of self-similar solutions : proof 
of Theorems 11.31 and 11.41 

Proof of Theorem ] 1.3\ : Let m > M and uq G Ym- 

From Theorem 15.11 we can fix e > such that is constant equal to L on 
[Af:,oo) with L < m. 

Pick to E (0, Tmax{ua)) ■ From Theorem l2.1l iv). we know that ^2,(^0, 0) > 0, then 

a — min — — > 0. Hence oq = min(l, — : ) > is small enough so 

^e(o.il ^ llK||oo,[o,i] 

that 



Let us set 



Ve{aox) < u{to,x) for all x E [0, 1] 
a{t) 



and 



(1 - eag qt)i 
ea^q 

Remark that a(0) — qq and 

d{t) =ea(i)i+'' forie [0,r*). 
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We shall show that 

u(t, x) — Ve{a{t)x) 

is a subsolution for {PDE)m with initial condition uo(^o)- 
Let us set T = mm{Tmaxiuo) - to,T*). 

Indeed, u{0,x) — Ve{aox) < uo{to,x) for all x £ [0, 1]. For all t G [0,T), we see 
that 0) = = uo(to + 1, 0) and u{t, 1) < L < m — uo(to + t,l). 
Moreover, a straightforward calculation shows that 



= 



where y = a{t)x. 

From the comparison principle (cf Lemma l2.ip . u(t) < uo{tQ+t) for alH £ [0,T). 
Now, if we assume that Tmaxiuo) > to + T* , then T — T* and by letting t go 
to T* , since a(t) — > +oo, we obtain 

L <uo{tQ + T*,x) for alla;e (0,1]. 

Since Uo{to + T* ,0) — 0, this contradicts the continuity of uo{to + T*) a,t x — 0. 
Hence, Tmax{uo) <ta + T*<oo. □ 

Proof of Theorem\1.4\ ■' From Theorem 15.11 we know the existence of ei > 
such that for e E (0, ei], < A + 1 and Ve is flat from x = A^. We set 
Ue,a(,(t,x) — Ve{a{t)x) whcrc a{t) = 3- and oq > A^. The calculation 

in the proof of the previous theorem and the fact that Vc is constant from 
the point x = A^ prove that (ue.ao)^ >^ is a family of solutions of problem 
(PDEm,) where M, = Ve(A + 1). 

Now, by using the same methods as in the proof of Lemma 15.61 we can prove 
that 

e i-> Ve is continuous from [0,ei] to C^([0,+oo)). (56) 

Then, e i->- = Ve{A + 1) is continuous in [0, ei] so its image is a compact 
interval /. Since Ue.ao blows up, we necessarily have > Af if e S (0, ei] so 
/ = [M, M+] with M+ > 0. 

Eventually, denoting = sup ^^ii£l^ n jg clear that Kf^ = sup ^"^^-^ since 

a;g(0,oo) x<£{a,A,] 

Ve is flat from x — A^. Then, we have 

AA[u,,„„(t)] = K^ait) - , , -T, 

qie-'iTmax - t) 1 

where 

1 



We now prove that is concave for small e, which implies that if^ = 1 so that 
the blow-up speed is known explicitly. Let e g [0, ei]. 

By Lemma there exists 6 > independent of e such that Ve > \ on [0,(5]. 

Since Ve is moreover nondecreasing on [0, 1], then min — > |. By (ISSl) . 

xe(a,A+i] ^ ^ 

there exists C > such that || KHoo < C* for all e G [0, ei]. 
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Let £2 = min(ei, ^ci-q ), e € (0, 62] and x G [0,A+1]. 
We have exVeix)'^'" < e2xC'^-i < fa; < Ve{x) and since 



x'^-'^Velx) = -V,{xy Ve{x) - exV,{xy-'i , 
we deduce that is concave and = Ve{0) = 1 for e e [0, €2]- 



□ 
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